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Abstract 

We apply the Thermal Field Theory methods to study the propagation of photons in a plasma 
layer, that is a plasma in which the electrons are confined to a two-dimensional plane sheet. 
We calculate the photon self-energy and determine the appropriate expression for the photon 
propagator in such a medium, from which the properties of the propagating modes are obtained. 
The formulas for the photon dispersion relations and polarization vectors are derived explicitly 
in some detail for some simple cases of the thermal distributions of the charged particle gas, and 
appropriate formulas that are applicable in more general situations are also given. 
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I. INTRODUCTION 



It is well known that when elementary particles, such as photons or neutrinos, propagate 
through a medium, the effects of the background particles influence their properties in 
important ways. Among the various approaches that exist to study the effects of a medium 
on the properties and interactions of particles, the methods of Thermal Field Theory (TFT) 
have proven to be very useful. Largely motivated by the original work of Weldon[l-3|, these 
methods have been applied to the problems mentioned above and other similar ones, and 
they have been helpful for understanding the physics involved and also from a computational 
point of view. 

The present work is concerned with similar calculations, but with the distinction that 
the medium consists a gas of particles that are confined to live in a two-dimensional plane 
sheet, or layer. A typical system of this type is an ordinary plasma in which the electrons 
are confined to a plane sheet, that we can take to be the z = plane. The quantities of 
interest are the usual ones such as, for example, the dispersion relations of the propagating 
photon modes, the damping and the transition rates for various processes. 

Before continuing we want to stress that this is not the same thing as what is usually called 
QEDs (or QED in 2+1 dimensions), which has been studied in the literature % l5|. QED3 
describes a system that, with regard to the space coordinates, has cylindrical symmetry 
and the physics, being independent of the z coordinate, can be studied by considering a 
two-dimensional cross section. Thus, for example, the electron in QED^ is really a line of 
charge in the three-dimensional world, and the Coulomb potential between two such electrons 
is logarithmic. In contrast, in the system we are considering, the electron is an ordinary 
point charge, which is confined to the z = plane, but the Coulomb potential between two 
electrons is the usual 1/r potential. 

The method that we apply here to study these systems could be useful in the context of 
astrophysical 6-8| as well as plasma physicsjsl and condensed matter 0, [ll| applications. 



and they are also interesting in their own right because they can be useful in the study of 
physical systems of current interest in which a plasma is confined to a layer [3] or a wire 13 1. 

In the context of TFT, the distinctive feature of the system that we are considering is that 
the medium is not isotropic over the (three-dimensional) space. Consequently, the thermal 
propagators that are used ordinarily in TFT calculations for the case of homogenous and 
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isotropic media, are not the appropriate ones for the present case. Therefore, in order to 
use the TFT methods to study the model of the two-dimensional plasma layer, a crucial 
requirement is finding the appropriate set of thermal propagators that must be used. 

In the present work we considered the simplest situation of an ordinary gas of electrons, 
which are confined to a plane sheet, but are otherwise free. Our main goal has been to 
formulate the TFT approach to the model of the two-dimensional plasma layer that we 
have described. The important steps to this end are taken in the first part of the present 
paper, where we determine the appropriate set of thermal propagators. The charged particle 
propagators are very similar in form to the standard three-dimensional form. But, as we will 
see, the propagation of the photon in the layer is described by en effective field which has a 
corresponding propagator that is very different from the usual one. The photon propagator 
is an important quantity because its inverse determines the bilinear part of the effective 
action or, equivalently, the equation of motion for the photon effective field, from which the 
dispersion relations and wave functions of the propagating photon modes can be obtained. 
As an application, in the second part we carry out a one-loop calculation of the photon 
self-energy in that medium and, as a specific example, we consider in detail the calculation 
of the longitudinal dispersion relation. There we compare our approach and results for this 
calculation with the results that are known for this system in literature jl^. IsL which have 
been obtained using the so-called static local field correction approximation[16]. There we 
show explicitly that our results for this calculation reduce to those known results when the 
appropriate limits are taken and/or approximations are made, which in particular involves 
approximating the photon propagator by its static (zero frequency) limit. However, the 
formulas that we obtain for the self-energy and dispersion relations in general, and the 
longitudinal dispersion in particular, can be used for a wider range of conditions in which 
those approximations (such as the static local field correction approximation) and limits are 
not justified. 

II. NOTATION AND KINEMATICS 

As usual, we denote by the velocity four-vector of the medium. Adopting the frame 
in which the medium is at rest, we set 

u'^ = {l,6), (2.1) 



and from now on all the vectors refer to that frame. Thus, denoting by n the unit vector 
perpendicular to the plane layer, we introduce the four-vector 

n^^{0,n), (2.2) 

and denote the momentum four-vector of a photon that propagates in the plane is by 

k'i^iuj,K^), (2.3) 



where 



It is useful to define 



where 



as well as the tensors 



K±-n^O. (2.4) 



u^, = g^.yu" , (2.5) 



Qf^u = Qnv p — , (2.6) 



It is useful to note that the tensors Rfjti, and Q^^i, are transverse to /c^, i.e., 

k'iR^, = k'iQ^, = Q, (2.8) 



and they satisfy the relations 



p pA p 

R^xQ\ = 0. (2.9) 



III. THE MODEL 

A. The electron free field 

We envisage a slab of area in the xy plane and thickness L^. The electrons are 
confined to live within that slab but are otherwise free to move within it. Eventually, the 
limits L — > oo and — > will be taken in a suitable way. 



Making use of the Furry picture, the electron field is expanded in terms of the one-particle 
wavefunctions which we take to be of the form 

H{z)e'P^-'^^ , (3.1) 

where f)± is the component of p in the xy plane and similarly for x^. The wavefunctions can 
be taken to satisfy periodic boundary conditions in the xy plane as usual. In principle the 
function H{z) is labeled by some quantum number. However, since we are contemplating 
taking the — )■ limit, for our purposes we assume that only the lowest lying state survives. 
Moreover, for the same reason, the particular form that it may take is not relevant, as long 
as it is consistent with our assumptions, and in particular we can adopt 

{1 for - ^ < 2; < ^ 
(3.2) 
otherwise 

Then, making the passage to the continuous momentum plane waves in the xy plane (L 



00), the electron field has the form 

ix^ = ( 



^(x) = ( ^] i,{x^) , (3.3) 



^{x^) = I [a(px, s)u{p^, s)e-'^--^- + h\p^, s)v{p^, s)e^^-^-] , (3.4) 



where 



with E = \jp1_ + m^, = (x°,x_l) and = {E,p±). The spinors u are the standard Dirac 
spinors normalized such that 

uu = 2m , (3.5) 
and the creation and annihilation operators satisfy 

{a{p^,s),a*{p'^,s')} = {2Tif2E5^^\p^-p'^)5s,s'^ (3.6) 

with analogous relations for the spinors v and the h operators. 

The passage to the planar layer will be made by taking the — > limit at the appropriate 
stage, which in essence reduces to use 

H{0) = 1, (3.7) 

where S{z) is the Dirac delta function. For practical purposes Eqs. f l3.3p -( l3l7|) can be taken 
as the equations that define our model for the plane layer. 



B. Electron thermal propagator 



From the point of view of TFT, the important ingredient for carrying out the calculation 
of the photon self-energy is the electron thermal propagator, which in turn can be determined 
by well established rules in terms of the free-field propagator. The propagator associated 
with the electron free-field ip± can be determined in various ways, in analogy with the usual 
case, and the results look similar. The various components of the thermal propagator matrix 
are simply 

1 



S±22{PA 
S±2liPA 



(|^_L + me) 
(;^_L + me) 



p'j_ — ml + ie 
-1 



+ 2ni6{pj_ - m^)rie{p± ■ u) 



p]_ — ml 



le 



+ 2Txi5{p\ - ml)r]e{p± ■ u) 



where 



with 



= + me)27ri [r7e(p± ■ u) - 9i-p± ■ u)] , 
= (|^_L + m^)2ni [r]eip± ■ u) - 9{p± ■ u)] , 

Veip) = 0{p ■ U)fe{p ■ U) + 0{-p ■ U)fe{-p " u) 



fe{x) 



(3.8) 



(3.9) 



g/3(2:-/Xe) _|_ I 

1 



(3.10) 



g/3e(a;+Aie) _|_ ]_ 

and 6{x) is the step function. Here /3e and /ie are the inverse temperature and the chemical 
potential of the electron gas, respectively. The total number of particles in the gas can be 
calculated from 



N 



(fx ip^^ifj 

(fp± 
(2^ 



Tr 



S±ii{p±)j^ 



(3.11) 



where we have set J dxdy — > L^. Using Eq. (13. 3 p and the formulas given above for the 
propagator, this yields 



N/L' 



rte + Ue 



where 



<fp± 



1 



(27r)2e/3(ET/.e) + i ' 

which represent the surface density of electrons and positrons respectively. 



(3.12) 
(3.13) 



IV. PHOTON PROPAGATION IN THE LAYER 



A. Photon effective field 

Using Eqs. f l3.3p and f l3.7p . the usual interaction Lagrangian term ■iJj'y^^iIjA^ yields the 
following term in the action 

S;^t = -^ j d'x^h^ijA,, (4.1) 

where 

K = ^mL=o • (4-2) 

This indicates that is the relevant electromagnetic field variable, and it is the one that we 
should focus on. Thus, we regard A^ as the effective field for the photon, and our goal is to 
determine its effective action, or equivalently its equation of motion, including the thermal 
corrections. 

Formally, what we want to do is to integrate out all the dynamical field variables except 
A^ itself. Following the usual functional method of quantization of the electromagnetic 
field [itI. and adapting it to the present model, a convenient way to proceed is to introduce 
in the action an external current 

J'{x) = , (4.3) 

with Jf^lx^) satisfying 

d±-J{x±) = 0, (4.4) 

where = (gf^, ~gf^)- Eq. fl4.4p ensures that we are selecting only the transverse (gauge 
invariant) part of A^, which is the physically meaningful one, since the longitudinal part 
decouples. The classical field, which we denote by A^^\ in the presence of both, the external 
current and the interaction given by S^^^^ in Eq. (14. ip . is then defined by 

4.5 
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B. Photon propagator 



Following the usual argument, and remembering Eq. (13. 7p . the generating functional for 
the free electromagnetic field in the layer is 

Z[J] oc exp j d^x±d^x'J^{x±)AFf,u{x± - x'^)r{x'^)^ , (4.6) 

where Ap^^{x± — x'j_) is obtained from the standard photon propagator Ap^^{x — x') by 
setting the coordinates normal to the plane {z and z') equal to zero. Therefore, taking into 
account Eq. (14. 4p . the propagator for the effective field in the plane is given, in momentum 
space, by 

Ap,.{kA_) = ~g,a~9up ( / ^ Af{k)^ , (4.7) 

where g^^ has been defined in Eq. (12.61) and, in the integrand, the momentum vector k is 
decomposed in the form 

kf, = k±f, + , (4.8) 
with k±^ as given in Eq. (12. 3p . Writing 

AFuuik) = „ ^^'^ — h gauge-dependent terms, (4.9) 
+ ie 

and carrying out the integration over Ky, we then obtain the propagator for the effective 
photon field as 

Ap^^uik^) = -A{k±) {R^.u + Qm-) , (4.10) 
where the tensors R and Q have been defined in Eq. (12.70 . and 

[ UJ>K 

A{k^) = (4.11) 

with 

k=\k±\. (4.12) 



C. Equation of motion 

Armed with the expression for the free photon propagator in the plane, the bilinear part 
of the effective action for Aj/^ is then given, in momentum space, by 

S^''> = / 1^ [D^'ik^) + ^'^(k^)] Aif\k^) - A('>{k^) ■ J{k^)] , (4.13) 
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FIG. 1: One-loop diagram for the photon thermal self-energy matrix, 
where D^'^{k±) is defined 

A^p\k^)Dx.ik^) = , (4.14) 

and TTfj^iy is the photon self-energy in the medium. The dispersion relations, and the cor- 
responding polarization vectors, of the propagating photon modes can be determined by 
solving the equation of motion for the classical field in the absence of the external current, 
that is 

[D^^''{k^) + nf^'^{k^)]A^^\k^) = 0. (4.15) 
D. Photon self-energy and dispersion relations 

We denote the components of the thermal self-energy matrix by tt^'^''^ , which to the lowest 
order are determined by calculating the one-loop diagram shown in Fig. [1] The physical 



self-energy function that appears in Eq. fl4.15l) is determined from the relations 2l| 



Re7r^.(fcx) = Re7r(i;)(A;^) 

Imfc^^k^) = (4.16) 

where 

= ^r- — r- (4-17) 
In general vr^,^ satisfies the transversality condition 

k^f^^^u = k'i7T^, = 0, (4.18) 



9 



as a consequence of the conservation of the electromagnetic current. Since vr^i/ is a function 
only of the vectors kj_ and u^, it follows that it can be expressed in the form 1| 

TT^^ = TItRi,,, + T^lQfiu ■ (4.19) 

In the more general case, n^^, can contain a term proportional toQ, 13] 

P^.u = -e^.,apu''ki. (4.20) 



However, this term is odd under the Parity and CP transformations and therefore it can 
arise only if the system is not invariant under those two transformations, such as in a gas of 
chiral molecules or if the weak interactions are included. Since we are not considering any 
of those effects here, for our purposes Eq. fl4.19p is the most general form. 
Using Eq. f l4.10p . V^^lkx) is given by 

D'^'^ik^) = -A-\k^) {R,, + Q^,) , (4.21) 

and, remembering Eq. f l4.19p . the dispersion relations 

A'\k^)-iTL{k^) = 0, (4.22) 
K-\k^) - TTT{k±) = (4.23) 

follow, where A{k±) is given in Eq. (14. lip . 



V. ONE-LOOP CALCULATION 

A. One- loop formula for the self- energy 

We want to apply the results of the one-loop calculation of the photon self-energy, to the 
dispersion relations given in Eqs. f l4.22p and fl4.23p . We will restrict ourselves to the real 
part of the dispersion relations, and therefore we consider here the calculation of nj^J^ , from 
which the real part of the physical self-energy is determined by means of Eq. (14.161) . 

Referring to Fig. [H 

m^^A;^) = (-l)(-0'Tr J ^^i,zSMp± + k±h.iSMP±) ■ (5-1) 

When the formula for S±n given in Eq. (13.80 is substituted in Eq. (15. ip . there are three types 
of terms. The term that contains two factors of r/e contributes only to the imaginary part 
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of the self-energy and, since we restrict ourselves here to the real part, we do not consider 
it further. The remaining terms then yield 



(5.2) 



where vr^^-* is the standard vacuum polarization term, while the background dependent con- 
tribution is given by 



(27r)22E 



ife + h) 



-'flU 



In this formula. 



p'i = {E,p^), E = ^Jpl+ml 

and /e,e denote the particle and antiparticle number density distributions given by 

1 



fe 



(5.3) 

(5.4) 
(5.5) 

(5.6) 



with the minus(plus) sign holding for the electrons(positrons), respectively. The integral in 
Eq. (15. 3p is to be interpreted in the sense of its principal value part. 

It is easily verified that vr^™) satisfies the transversality condition Eq. (I4.18p . and therefore 
it can be decomposed in the form 



TT 



(m) 



as already indicated. The functions 7r!f^l can be found by projecting Eq. (15. Sp 
tensors R^u and Qfj,u, and this procedure yields 



(5.7) 
with the 



(m) 



TT 



(m) 



4e^^5e, 



(5.^ 



where 

A, 
5, 



d^Pi. 

(Ppi. 
(27r)22E 



{fe + h) 
ife + h) 



2ml ~ 2p_L ■ fcj 
k\ + 2p^ ■ k^ 



+ {k^ -kx 



2{p± -uy + 2{p±-u){k±-u) -p±-k± 
kl + 2p^ ■ k^ 



+ {k± -kj 



(5.9) 
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These two integrals are very similar to the corresponding ones that appear in the three- 
dimensional case, which have been analyzed in the literature in considerable detail, for 
various limiting values of the photon momentum and several conditions of the charged 
particle gas. The methods employed there are applicable here as well. While their evaluation 
is not possible for the general form of the distribution functions, some useful results can be 
obtained by considering special cases. 



B. Low momentum limit 



In many situations of interest, the photon momentum is such that 

w, K<£'e, (5.10) 
where E^, is the typical energy of the electrons in the gas. In this case, borrowing the results 



= 5e + ^ /tTT^ ^^ h^(/e + /.), (5.11) 



given in Eq. (3.2) of Ref. |20|, 

_1 r fpj^ f Vp-K \ _d_ 
^ ~ 2] (27r)2 \u-v,-^) dE^^'^^~'' 
u r (fp± / \ d 

2 J (27r)2 \uj-Vp-h) dE' 

with Vp = p±/E denoting the velocity of the particles in the background. The above formula 

for Be can be rewritten by multiplying the integrand by the factor 

— {oj — Vp ■ K + Vp ■ K.) . (5-12) 

CO 

The first two terms integrate to zero, while the third one leads to 

Be = --I^A S^eSA _^(/^ + /,) . (5.13) 

We stress that the expressions given in Eq. fl5.1ip are derived from Eq. fl5.9p by expanding 
the integrands in terms oik/E and retaining only the terms that are dominant when k/E ^ 
0. For a non-relativistic gas, Eq. (15.1 ip holds for u, k, ^ rrie. If the gas is extremely 
relativistic, Eq. (15. lip also holds for u, k > nie, subject to Eq. (I5.10p . 

Furthermore, up to this point no assumption has been made regarding the nature of the 
electron background. Apart from the restriction on k, Eqs. (15. lip and (I5.13P hold for a 
relativistic or non-relativistic gas, whether it is degenerate or not. Accordingly, they serve 
as a convenient starting point to find the dispersion relations provided that Eq. (I5.10p is 
verified. As an example we consider one specific situation below. 
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C. Longitudinal dispersion relation 



We consider the case of a non-relativistic electron gas, in the regime in which Eq. f lS.lOp 
is vahd, and seek the real solutions to the longitudinal dispersion relations in the long 
wavelength limit, 

Cc;>Kfe, (5.14) 

where Ve is the typical velocity of the electrons. Then, neglecting f-g in Eq. fl5.13p . the 
calculation of B^. in this limit gives 

and using Eqs. fl4J6|) . (Q and flS^Sjl . 



Be = ^. (5.15) 



Revr^ = i-^p-«:2) . (5.16) 
A real solution of Eq. (14.221) is obtained for a; < k by solving 



which gives 

.2\l/2 



2 -^2 



where 



(5.1^ 



a = -. (5.19) 

2m. ^ ' 



The function is depicted in Fig. |2l It follows from Eq. (I5.18P that ~ k for k ^ a, 
while for k ^ a it is given approximately by 

u'^ = aK. (5.20) 

This solution is valid as long as 

«:<4' (5-21) 

so that the condition u ^ nve is satisfied and therefore Eq. (15.151) remains valid. 

This result, which implies that in this example system the plasma frequency is 
momentum-dependent, is known in the plasma physics and condensed matter literature jl^. 



15|. In those contexts, it is derived using the static local field correction approximation, 
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FIG. 2: Plot of the longitudinal dispersion relation given in Eq. ()5.18p . 



in which the dynamic corrections are replaced by their static [u — )■ 0) values IGj. In our 
notation, that approximation corresponds to using the expressions 

A{k^) ~ lim A(/cx) = 7^ , 

RevTL ~ Rem = 5.22 

in Eq. f l4.22p . The solution of the equation thus obtained is indeed given by Eq. f l5.20p . But, 
as we have seen, Eq. (15.200 is an approximation to the solution given in Eq. fl5.18p . which 
we have obtained by using the complete expressions for the propagator and self-energy. 

Therefore, in our case, by using the complete inverse propagator in Eq. (14.220 . we are 
able to obtain the dispersion relation given in Eq. (I5.18p . which is valid for all values of k, 
including k ^ a which are outside the range of validity (k ^ a) of the static local field 
correction approximation. 

Our motivation for going through this exercise here is to evidence that some known results 
are not only reproduced, but also generalized, using the method we plan to use, which in 
turn indicate that they can be applied systematically to study the problems that we have 
mentioned in this type of system. 
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VI. CONCLUSIONS 



In the present work we have formulated the TFT approach to the model of the two- 
dimensional plasma layer, that is, a system in which the electrons are confined to a plane 
sheet. As emphasized in the Introduction, this is not equivalent to what is usually called 
QED^ (or QED in 2+1 dimensions) which describes a two-dimensional cross section of a 
system that has cylindrical symmetry and the physics is independent of the z coordinate. 
Thus, for example, while the electrons in QED3 are really lines of charge in the three- 
dimensional world and the Coulomb potential between them is logarithmic, in the system 
we are considering the electrons are ordinary point charges, which are confined to the z = 
plane, but the Coulomb potential between them is the ordinary 1/r potential. 

An important step in this direction was to determine the appropriate set of thermal 
propagators. As an application, we performed the one-loop calculation of the photon self- 
energy in that medium, and we considered in detail the calculation of the longitudinal photon 
dispersion relation. We made contact with previous calculations of that quantity, that had 
been obtained using other approaches, and in particular we showed that our results reduce to 
those known results when the appropriate limits are taken and/or approximations are made. 
However, the formulas we obtained are more general and can be used for a wider range of 
conditions in which those approximations and limits are not justified. We considered the 
simplest situation of an ordinary gas of electrons, which are confined to a plane sheet, but 
are otherwise free. However the method allows us to consider variations of the model in a 
systematic way, such as the effects of anisotropics and/or external fields. 

There is a considerable amount of literature on the applications of Thermal Field Theory 
(TFT) methods to the study of the properties of the electron gas (plasmon dispersion rela- 
tions and so on) in the three-dimensional space. Although those systems have been studied 
by a variety of methods, the application of TFT methods to study them have been useful 
in many ways. 

We believe that the applications to the two-dimensional sheet that we have described, are 
novel TFT calculations which are expected to be equally attractive and useful. For example, 
they can have astrophysical applicationsjo-Sl; they establish a good point of contact between 
the applications of TFT methods and dominant themes in plasma physics (such as the effects 
of anisotropics and instability studies); they may be are relevant for laboratory experiments 
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lll. 



in fundamental disciplines of physics such as plasma physics jol and condensed matter 

Such calculations are also interesting in their own right because they can be useful in the 
study of p hysical systems of current interest in which a plasma is confined to a layer 12| or 



a wire 



13j, and they may also be helpful for studying and understanding analogous effects 



and issues in more complicated systems such as the non-abelian plasmas. 
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